Simplified S1 vertebral bone quality score independently predicts proximal junctional kyphosis after surgery for degenerative lumbar scoliosis.

OBJECTIVE: Our study aimed to assess the effectiveness of the simplified S1 vertebral bone quality (VBQ) score in predicting the incidence of proximal junctional kyphosis (PJK) after surgery for degenerative lumbar scoliosis (DLS). METHODS: We reviewed 122 patients with DLS who underwent posterior lumbar decompression and long-segment fusion surgery in our hospital from January 2016 to December 2020. The patients were classified into PJK group and non-PJK group. S1 VBQ scores are determined by signal intensity measurements taken from the mid-sagittal plane of T1-weighted non-contrast MRI. Logistic regression analysis was used to identify factors associated with PJK. Receiver-operating characteristic curve (ROC) analysis was used to evaluate the value of S1 VBQ score in predicting pedicle PJK after DLS. RESULTS: 122 DLS patients (90 females and 32 males) met the inclusion criteria. In addition, 27 patients (22.13%) had PJK at the time of last follow-up. VBQ was higher in PJK patients than non-PJK patients (3.58 ± 0.67 vs. 3.08 ± 0.54, p < 0.001). Preoperatively, patients in the PJK group had a greater TLK than those in the non-PJK group (20.00 ± 6.22 vs. 16.86 ± 5.38, p = 0.011). After surgery, patients in the PJK group had greater TLK (p < 0.001) and PJA (p < 0.001) compared with the non-PJK group. At final FU, patients in the PJK group had greater TK (p = 0.002), TLK (p < 0.001), SVA (p < 0.001), and PJA (p < 0.001) than patients in the non-PJK group (Table 4). In multivariate logistic regression analysis, higher VBQ score (OR 4.565, 95% CI 1.43-14.568, p = 0.010), advanced age (OR 1.119, 95% CI 1.021-1.227, p = 0.016), and larger TLK (OR 1.191, 95% CI 1.041-1.362, p = 0.011) were significant predictors of postoperative PJK in patients with DLS (Table 6). A statistically significant positive correlation existed between VBQ score and PJA change (r = 0.370, p < 0.001). We created ROC curves for VBQ scores as predictors of PJK with a diagnostic accuracy of 72.1% (95% CI 60.15-82.9%.The ideal limit for the VBQ score was 3.205 (sensitivity: 77.8%, specificity: 81.4%). CONCLUSION: To the best of our knowledge, this is the first study to evaluate the effectiveness of the S1 VBQ score in predicting postoperative PJK in DLS. Our study included major risk factors and found that S1 VBQ score was a significant predictor of PJK in patients undergoing DLS surgery. The higher the S1 VBQ score, the higher the probability of PJK.

Cryptic diversity begets challenges and opportunities in biodiversity research.

How many species of life are there on Earth? This is a question that we want to know but cannot yet answer. Some scholars speculate that the number of species may reach 2.2 billion when considering cryptic diversity and that each morphology-based insect species may contain an average of 3.1 cryptic species. With nearly two million described species, such high estimates of cryptic diversity would suggest that cryptic species are widespread. The development of molecular species delimitation has led to the discovery of a large number of cryptic species, and cryptic biodiversity has gradually entered our field of vision and attracted more attention. This paper introduces the concept of cryptic species, how they evolve, and methods by which they may be discovered and confirmed, and provides theoretical and methodological guidance for the study of hidden species. A workflow of how to confirm cryptic species is provided. In addition, the importance and reliability of multi-evidence-based integrated taxonomy are reaffirmed as a way to better standardize decision-making processes. Special focus on cryptic diversity and increased funding for taxonomy is needed to ensure that cryptic species in hyperdiverse groups are discoverable and described. An increased focus on cryptic species in the future will naturally arise as more difficult groups are studied, and thereby, we may finally better understand the rules governing the evolution and maintenance of cryptic biodiversity.

Electronic structure engineering on Co-based metal-organic frameworks for concurrent electrocatalytic hydrogen generation and formate electrosynthesis.

A facile one-step solvothermal method was developed to synthesize Ir-doped Co-based metal-organic framework (CoIr-MOF) nanoarrays as a bifunctional electrocatalyst for water-glucose co-electrolysis. It was demonstrated that in situ incorporation of a low-content of Ir cations could modulate the electronic structure of Co active centers and thus boost the electrocatalytic performance towards both the hydrogen evolution reaction and glucose-to-formate oxidation reaction.

Comparison of nerve versus tendon transfer for radial nerve palsy.

OBJECTIVE: This study aims to investigate the choice of intervention time and operation mode between nerve and tendon transfer for the treatment of radial nerve palsy (RNP). METHODS: 46 RNP patients underwent nerve transfer (n = 22) and tendon transfer (n = 24). The intraoperative blood loss, main incision length, operation duration, and length of hospital stay and follow-up period of patients in these two groups were recorded and compared. The range of motion (ROM) of the elbow, wrist, fingers, and thumb, the hand grip and pinch strength, the Disabilities of Arm, Shoulder, and Hand (DASH) and the 36-Item Short Form Health Survey (SF-36) scores were measured and compared preoperatively and postoperatively between the two groups. RESULTS: The ROM of thumb and the hand grip strength of patients in the nerve transfer group were greater than that in the tendon transfer (P < 0.05). Both of the two groups indicated significant improvements in the ROM of elbow, wrist, finger, thumb and the hand grip and pinch strength (P < 0.05) postoperatively. The DASH scores decreased significantly at 6 months (P < 0.05) and 12 months (P < 0.05) after surgery in both groups, while the postoperative SF-36 scores significantly increased (P < 0.05). There was no significant difference in postoperative DASH and SF-36 scores between the two groups (P > 0.05). CONCLUSION: In summary, both nerve and tendon transfer techniques are effective treatments for RNP. Nerve transfer is particularly advantageous for early RNP, while tendon transfer is suitable for patients with radial nerve injury more than one year.

Tree communities and functional traits determine herbivore compositional turnover.

There are many factors known to drive species turnover, although the mechanisms by which these operate are less clear. Based on comprehensive datasets from the largest tree diversity experiment worldwide (BEF-China), we used shared herbivore species (zeta diversity) and multi-site generalized dissimilarity modelling to investigate the patterns and determinants of species turnover of Lepidoptera herbivores among study plots across a gradient in tree species richness. We found that zeta diversity declined sharply with an increasing number of study plots, with complete changes in caterpillar species composition observed even at the fine spatial scale of our study. Plant community characteristics rather than abiotic factors were found to play key roles in driving caterpillar compositional turnover, although these effects varied with an increasing number of study plots considered, due to the varying contributions of rare and common species to compositional turnover. Our study reveals details of the impact of phylogeny- and trait-mediated processes of trees on herbivore compositional turnover, which has implications for forest management and conservation and shows potential avenues for maintenance of heterogeneity in herbivore communities.

The impact of heavy metal pollution on wild bee communities in smallholder farmlands.

Wild bees provide important pollination services, but they face numerous stressors that threaten them and their ecosystem services. Wild bees can be exposed to heavy metal pollution through the consumption of nectar, pollen, and water, which might cause bee decline. While some studies have measured heavy metal concentrations in honeybees, few studies have monitored heavy metal concentrations in wild bees or explored their potential effects on wild bee communities. To investigate the impact of heavy metal pollution on wild bee communities, heavy metal concentrations, including vanadium (V), chromium (Cr), nickel (Ni), cadmium (Cd), Zinc (Zn) and lead (Pb) in multiple wild bee species were measured. Multiple wild bee species, including: Xylocopa tranquabaroroum, Eucera floralia, Apis cerana, and small bee mixtures (representing multiple small wild bee species) were sampled from 18 sites in Quzhou, Zhejiang Province, China. The findings demonstrated that there were significant differences in heavy metal concentrations among different bee species. The concentrations of V, Zn, Cd, and Pb in X. tranquabaroroum, the largest bee species in this study, were lower than that in the other three sample groups. Furthermore, there were significant negative correlations between heavy metal pollution and wild bee diversity and species richness, but not with abundance. Particularly, there was no significant relationship between heavy metal pollution and the abundance of small bees. Given these worrying findings, monitoring multiple heavy metals in wild bees should be conducted for protecting wild bee diversity and securing their pollination services.

Phylogeny-based assignment of functional traits to DNA barcodes outperforms distance-based, in a comparison of approaches.

The full potential for using DNA barcodes for profiling functional trait diversity has yet to be determined in plants and animals; thus, we outline a general framework for quantifying functional trait diversity of insect community DNA and propose and assess the accuracy of three methods for achieving this. We built a novel dataset of traits and DNA barcodes for wild bees in China. An informatics framework was developed for phylogeny-based integration of these data and prediction of traits for any subject barcodes, which was compared with two distance-based methods. For Phylogenetic Assignment, we additionally conducted a species-level analysis of publically available bee trait data. Under the specimen-level dataset, the rate of trait assignment was negatively correlated with distance between the query and the nearest trait-known reference, for all methods. Phylogenetic Assignment was found to perform best under several criteria; particularly, it had the lowest false-positive rate (rarely returning a state prediction where success was unlikely; where the distance from query to the nearest reference was high). For a wider range of compiled traits, conservative life-history traits showed the highest rates of assignment; for example, sociality was predicted with confidence at 53%, parasitism at 44% and nest location at 33%. As outlined herein, automated trait assignment might be applied at scale to either barcodes or metabarcodes. With further compilation and databasing of DNA barcode and trait data, the rate and accuracy of trait assignment is expected to increase to the point of being a widely viable and informative approach.

A Chromosome-level Genome Assembly and Evolution Analysis of Andrena camellia (Hymenoptera: Andrenidae).

Andrena camellia, an effective pollinator of the economicallysignificant crop Camellia oleifera, can withstand the toxic pollen of C. oleifera, making A. camellia a crucial for resource conservation and cultivation of C. oleifera. In this study, the whole genome of A. camellia was sequenced on the Oxford Nanopore platform. The assembled genome size was 340.73 Mb including 50 scaffolds (N50=47.435 Mb) and 131 contigs (N50=17.2 Mb). A total of 11, 258 protein-coding genes were annotated, in addition, 1,104 non-coding RNAs were identified. Further analysis that some chromosomes of A. camellia have a high level of synteny with those of Apis mellifera, Osmia bicornis and Andrena minutula. Thus, our reported genome of A. camellia serves as a valuable resource for studying species evolution, behavioral biology, and adaption to toxic pollen of C. oleifera.

Yttrium-doped cobalt-based metal-organic frameworks for selective electrooxidation of glycerol coupled with hydrogen production.

A facile strategy was developed for one-step synthesis of yttrium(Y)-doped cobalt-based metal-organic framework (CoY-MOF) nanosheet arrays. It was demonstrated that in situ doping with a low content of Y can tailor the electronic structure of the MOF structures, thereby improving their electrocatalytic performance for both hydrogen evolution and glycerol oxidation reactions.

K-line tilt as a novel potential risk factor for cervical Modic change: a retrospective study.

BACKGROUND: Cervical sagittal parameters are important parameters that reflect the mechanical stress in the sagittal plane of the cervical spine and are an important basis for predicting the clinical status and prognosis of patients. Although it has been confirmed that there is a significant correlation between cervical Modic changes and some sagittal parameters. However, as a newly discovered sagittal parameter, there is no report on the relationship between the K-line tilt and the Modic changes of cervical spine. METHODS: A retrospective analysis was performed for 240 patients who underwent cervical magnetic resonance imaging scan for neck and shoulder pain. Among them, 120 patients with Modic changes, namely the MC(+) group, were evenly divided into three subgroups of 40 patients in each group according to different subtypes, namely MCI subgroup, MCII subgroup and MCIII subgroup. One hundred twenty patients without Modic changes were included in MC(-) group. We measured and compared the sagittal parameters of cervical spine among different groups, including K-line tilt, C2-C7 sagittal axial vertical distance (C2-C7 SVA), T1 slope and C2-7 lordosis. Logistic regression was used to analyse the risk factors of cervical Modic changes. RESULTS: The K-line tilt and C2-7 lordosis were significantly different between MC(+) group and MC(-) group (P < 0.05). The K-line tilt greater than 6.72° is a risk factor for Modic changes in cervical spine (P < 0.05). At the same time, the receiver operating characteristic curve showed that this change had moderate diagnostic value when the area under the curve was 0.77. CONCLUSION: This study shows that the K-line tilt greater than 6.72° is a potential risk factor for Modic changes in cervical spine. When the K-line tilt is greater than 6.72°, we should be alert to the occurrence of Modic changes. TRIAL REGISTRATION NUMBER: 2022ER023-1.

The Whole Genome of Colletes collaris (Hymenoptera: Colletidae): An Important Step in Comparative Genomics of Cellophane Bees.

Gaining knowledge on bees is of the utmost importance due to the paramount role that they play in angiosperm pollination. Herein, we provide the first genome assembly of Colletes collaris, a pan-Eurasian cellophane bee. We sequenced 50.53 Gbp of long-read data plus 57.36 Gbp of short-read data in Oxford Nanopore Technologies and Illumina platforms, respectively. The genome assembly consisted of 374.75 Mbp distributed across 374 contigs, with L50 and N50 of 9 and 8.96 Mbp, respectively. We predicted the genome to comprise 20,399 protein-coding genes, 467,947 repeats, and 4,315 non-coding RNA genes. The transcriptome and mitochondrial genome of the species were also assembled. Gene family analysis with 15 insect species identified 14,417 families, 9,517 of them found in C. collaris. A dated phylogenomic analysis revealed high numbers of orthogroups experiencing rapid evolution within Colletes.

Impacts of Taxon-Sampling Schemes on Bayesian Tip Dating Under the Fossilized Birth-Death Process.

Evolutionary timescales can be inferred by molecular-clock analyses of genetic data and fossil evidence. Bayesian phylogenetic methods such as tip dating provide a powerful framework for inferring evolutionary timescales, but the most widely used priors for tree topologies and node times often assume that present-day taxa have been sampled randomly or exhaustively. In practice, taxon sampling is often carried out so as to include representatives of major lineages, such as orders or families. We examined the impacts of different densities of diversified sampling on Bayesian tip dating on unresolved fossilized birth-death (FBD) trees, in which fossil taxa are topologically constrained but their exact placements are averaged out. We used synthetic data generated by simulations of nucleotide sequence evolution, fossil occurrences, and diversified taxon sampling. Our analyses under the diversified-sampling FBD process show that increasing taxon-sampling density does not necessarily improve divergence-time estimates. However, when informative priors were specified for the root age or when tree topologies were fixed to those used for simulation, the performance of tip dating on unresolved FBD trees maintains its accuracy and precision or improves with taxon-sampling density. By exploring three situations in which models are mismatched, we find that including all relevant fossils, without pruning off those that are incompatible with the diversified-sampling FBD process, can lead to underestimation of divergence times. Our reanalysis of a eutherian mammal data set confirms some of the findings from our simulation study, and reveals the complexity of diversified taxon sampling in phylogenomic data sets. In highlighting the interplay of taxon-sampling density and other factors, the results of our study have practical implications for using Bayesian tip dating to infer evolutionary timescales across the Tree of Life. [Bayesian tip dating; eutherian mammals; fossilized birth-death process; phylogenomics; taxon sampling.].

Wideband Waveform Design for Distributed Precision Jamming.

Precision electronic warfare is a hot direction for future jamming technology development, and distributed precision jamming (DIPJ) is one of its typical application scenarios. The task objective of DIPJ is to design jamming waveforms so that the jamming energy generated by a set of ultra-sparse array transmitters can be focused in the jamming region of interest while being suppressed in other specific protected regions, which can be viewed as a distributed multiple-input and multiple-output system waveform design problem under a three-dimensional scenario. This paper extends the jamming signal model in DIPJ from narrowband to wideband based on previous work to address a broader range of jamming tasks. After extending the model to wideband signals, a method based on the traditional maximum total energy difference optimization objective is first given for comparison. A wideband jamming waveform design method based on the majorization minimization algorithm with the desired power spectrum matching as the optimization target is designed for the problem that the maximum energy difference method cannot focus energy well in the jamming region. The simulation results show that the presented method can make the jamming energy well concentrated in the target region and evenly distributed over the whole bandwidth, while the energy in the whole bandwidth is suppressed in the protected region.

The genome of sheep ked (Melophagus ovinus) reveals potential mechanisms underlying reproduction and narrower ecological niches.

BACKGROUND: Melophagus ovinus is considered to be of great veterinary health significance. However, little is known about the information on genetic mechanisms of the specific biological characteristics and novel methods for controlling M. ovinus. RESULTS: In total, the de novo genome assembly of M. ovinus was 188.421 Mb in size (330 scaffolds, N50 Length: 10.666 Mb), with a mean GC content of 27.74%. A total of 13,372 protein-coding genes were functionally annotated. Phylogenetic analysis indicated that the diversification of M. ovinus and Glossina fuscipes took place 72.76 Mya within the Late Cretaceous. Gene family expansion and contraction analysis revealed that M. ovinus has 65 rapidly-evolving families (26 expansion and 39 contractions) mainly involved DNA metabolic activity, transposases activity, odorant receptor 59a/67d-like, IMD domain-containing protein, and cuticle protein, etc. The universal and tightly conserved list of milk protein orthologues has been assembled from the genome of M. ovinus. Contractions and losses of sensory receptors and vision-associated Rhodopsin genes were significant in M. ovinus, which indicate that the M. ovinus has narrower ecological niches. CONCLUSIONS: We sequenced, assembled, and annotated the whole genome sequence of M. ovinus, and launches into the preliminary genetic mechanisms analysis of the adaptive evolution characteristics of M. ovinus. These resources will provide insights to understand the biological underpinnings of this parasite and the disease control strategies.

Functional and phylogenetic relationships link predators to plant diversity via trophic and non-trophic pathways.

Human-induced biodiversity loss negatively affects ecosystem function, but the interactive effects of biodiversity change across trophic levels remain insufficiently understood. We sampled arboreal spiders and lepidopteran larvae across seasons in 2 years in a subtropical tree diversity experiment, and then disentangled the links between tree diversity and arthropod predator diversity by deconstructing the pathways among multiple components of diversity (taxonomic, phylogenetic and functional) with structural equation models. We found that herbivores were major mediators of plant species richness effects on abundance, species richness, functional and phylogenetic diversity of predators, while phylogenetic, functional and structural diversity of trees were also important mediators of this process. However, the strength and direction differed between functional, structural and phylogenetic diversity effects, indicating different underlying mechanisms for predator community assembly. Abundance and multiple diversity components of predators were consistently affected by tree functional diversity, indicating that the variation in structure and environment caused by plant functional composition might play key roles in predator community assembly. Our study highlights the importance of an integrated approach based on multiple biodiversity components in understanding the consequences of biodiversity loss in multitrophic communities.

An integrative bioinformatics pipeline shows that honeybee-associated microbiomes are driven primarily by pollen composition.

The microbiomes associated with bee nests influence colony health through various mechanisms, although it is not yet clear how honeybee congeners differ in microbiome assembly processes, in particular the degrees to which floral visitations and the environment contribute to different aspects of diversity. We used DNA metabarcoding to sequence bacterial 16S rRNA from honey and stored pollen from nests of 4 honeybee species (Apis cerana, A. dorsata, A. florea, and A. laboriosa) sampled throughout Yunnan, China, a global biodiversity hotspot. We developed a computational pipeline integrating multiple databases for quantifying key facets of diversity, including compositional, taxonomic, phylogenetic, and functional ones. Further, we assessed candidate drivers of observed microbiome dissimilarity, particularly differences in floral visitations, habitat disturbance, and other key environmental variables. Analyses revealed that microbiome alpha diversity was broadly equivalent across the study sites and between bee species, apart from functional diversity which was very low in nests of the reclusive A. laboriosa. Turnover in microbiome composition across Yunnan was driven predominantly by pollen composition. Human disturbance negatively impacted both compositional and phylogenetic alpha diversity of nest microbiomes, but did not correlate with microbial turnover. We herein make progress in understanding microbiome diversity associated with key pollinators in a biodiversity hotspot, and provide a model for the use of a comprehensive informatics framework in assessing pattern and drivers of diversity, which enables the inclusion of explanatory variables both subtly and fundamentally different and enables elucidation of emergent or unexpected drivers.

Chromosome-Level Genome Assembly of Anthidium xuezhongi Niu & Zhu, 2020 (Hymenoptera: Apoidea: Megachilidae: Anthidiini).

Anthidiini, a large bee tribe characterized by light-colored maculations, represents nearly 1,000 pollinator species, but no genomes are yet available for this tribe. Here, we report a chromosome-level genome assembly of Anthidium xuezhongi collected from the Tibetan Plateau. Using PacBio long reads and Hi-C data, we assembled a genome of 189.14 Mb with 99.94% of the assembly located in 16 chromosomes. Our assembly contains 23 scaffolds, with the scaffold N50 length of 12.53 Mb, and BUSCO completeness of 98.70% (n = 1,367). We masked 25.98 Mb (13.74%) of the assembly as repetitive elements, identified 385 noncoding RNAs, and predicted 10,820 protein-coding genes (99.20% BUSCO completeness). Gene family evolution analyses identified 9,251 gene families, of which 31 gene families experienced rapid evolution. Interspecific chromosomal variation among A. xuezhongi, Bombus terrestris, and Apis mellifera showed strong chromosomal syntenic relationships. This high-quality genome assembly is a valuable resource for evolutionary and comparative genomic analyses of bees.

[Study on the safety and effectiveness of low-dose tranexamic acid in operation of multi-level continuous thoracic ossification of ligament flavum].

OBJECTIVE: To investigate the safety and effectiveness of low-dose tranexamic acid (TXA) in operation of multi-level continuous thoracic ossification of ligament flavum (TOLF). METHODS: A clinical data of 26 patients who underwent operation for multi-level continuous TOLF and met the selection criteria between July 2015 and January 2019 was retrospectively analyzed. Among them, 13 cases (group A) were received intravenous infusion of TXA (10 mg/kg) at 15 minutes before operation, and maintained the infusion at 1 mg/(kg·h) until the end of the operation; 13 cases (group B) were received the same dose of normal saline before and during operation. There was no significant difference in gender, age, body mass index, diseased segment, and preoperative hemoglobin, platelet count, activated partial thromboplastin time, prothrombin time, international normalized ratio (INR) between the two groups ( P>0.05). The hemoglobin, platelet count, activated partial thromboplastin time, prothrombin time, INR, the number of deep vein thrombosis of the lower extremities, operation time, intraoperative blood loss, postoperative drainage volume, total blood loss, and the time of drainage tube extubation in the two groups were recorded and compared. RESULTS: All operations in the two groups were successfully completed. Compared with group B, the operation time and time of drainage tube extubation in group A were shortened, and the intraoperative blood loss, postoperative drainage volume, and total blood loss were reduced. The differences between the two groups were significant ( P<0.05). None of the two groups received blood transfusion, and the hemoglobin level of group A at 24 hours after operation was significantly higher than that of group B ( t=5.062, P=0.000). The incisions in both groups healed and sutures were removed within 2 weeks after operation, and no complications occurred. There was no significant difference between the two groups in activated partial thromboplastin time, prothrombin time, INR, and platelet count at 24 hours after operation ( P>0.05). CONCLUSION: In multi-level continuous TOLF operation, intravenous administration of low-dose TXA can effectively reduce blood loss, shorten postoperative drainage time, and does not increase the risk of complications.

The complete mitochondrial genome of Cerceris quinquefasciata (Hymenoptera: Crabronidae).

The complete mitochondrial genome of the Cerceris quinquefasciata (Rossi, 1792) (Hymenoptera: Crabronidae) was obtained via next-generation sequencing. This mitochondrial genome is 16,188 bp in length with 37 classical eukaryotic mitochondrial genes and two A + T-rich region. All the 13 PCGs begin with typical ATN codons. Among them, eleven PCG genes terminate with TAA, two with T-. All of the 22 tRNA genes, ranging from 58 to 72 bp with typical cloverleaf structure except for trnS1, whose dihydrouridine (DHU) arm forms a simple loop. Phylogenetic analysis highly supported Crabronidae is the sister group of anthophila bees.

Tree diversity and functional leaf traits drive herbivore-associated microbiomes in subtropical China.

Herbivorous insects acquire microorganisms from host plants or soil, but it remains unclear how the diversity and functional composition of host plants contribute to structuring herbivore microbiomes. Within a controlled tree diversity setting, we used DNA metabarcoding of 16S rRNA to assess the contribution of Lepidoptera species and their local environment (particularly, tree diversity, host tree species, and leaf traits) to the composition of associated bacterial communities. In total, we obtained 7,909 bacterial OTUs from 634 caterpillar individuals comprising 146 species. Tree diversity was found to drive the diversity of caterpillar-associated bacteria both directly and indirectly via effects on caterpillar communities, and tree diversity was a stronger predictor of bacterial diversity than diversity of caterpillars. Leaf toughness and dry matter content were important traits of the host plant determining bacterial species composition, while leaf calcium and potassium concentration influenced bacterial richness. Our study reveals previously unknown linkages between trees and their characteristics, herbivore insects, and their associated microbes, which contributes to developing a more nuanced understanding of functional dependencies between herbivores and their environment, and has implications for the consequences of plant diversity loss for trophic interactions.